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ABSTRACT 

Aims. To study the formation of calcium carbonate, through the solid-gas interaction of amorphous Ca-silicate with 
gaseous C02, at elevated pressures, and link this to the possible presence of calcium carbonate in a number of circum- 
stellar and planetary environments. 

Methods. We use in-situ synchrotron X-Ray powder diffraction to obtain detailed structural data pertaining to the 
formation of the crystalline calcium carbonate phase vaterite and its evolution with temperature. 

Results. We found that the metastable calcium carbonate phase vaterite was formed alongside calcite, at elevated CO2 
pressure, at room temperature and subsequently remained stable over a large range of temperature and pressure. 
Conclusions. We report the formation of the calcium carbonate mineral vaterite whilst attempting to simulate carbon- 
ate dust grain formation in astrophysical environments. This suggests that vaterite could be a mineral component of 
carbonate dust and also presents a possible method of formation for vaterite and its polymorphs on planetary surfaces. 
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1. Introduction 

1.1. Carbonates in astrophysical environments 

On Earth, carbonates are the main constituent in chem- 
ically precipitated sedimentary rocks (e.g. limestone, 
dolomite) and are produced organically by marine organ- 
isms whose shells or skeletal structures are built out of cal- 
cium carbonate; however, the inorganic formation of car- 
bonates in terrestrial environments has been found to occur 
exclusively through the interaction of silicate minerals with 
liquid water containing dissolved COl~ . However, features 
indicative of carbonate minerals have been detected in a 
wide variety of astrophysical envi ronments, from planet ary 
bodies within the Solar System (Ehlmann et al. 2008 ) to 
dust particles surrounding stellar objects ( |Kemper et al. 
' ToppMli||2005| piiavassa et a"L]|2005[ ). 



2002 



1.1.1. Circumstellar Environments 

Carbonates have been detected in a number of circumstellar 
environments, exhibiting far-infrared (FIR) features around 
62/im and 92/im. The 92/im feature has been identified as 
calcite (CaCOs), being the only dust species known to have 
a strong feature at this wavelength. Dolomite (CaMgCOa) 
has a band at ~62/j,m, and has been suggested as a possible 
dust species; however due to other strong features of II2O 
and diopside (CaMgSi206) between 62-65/im this identi- 
fication is tenuous. Mg-rich and Fe-rich carbonates such 
as magnesite and siderite do not exhibit strong features in 
the far infrared and therefore cannot be confirmed as dust 
components. 



Kemper et al. (2002) report the detection of carbon- 



NGC6537. They argue that the carbonates in these envi- 
ronments require formation by non-aqueous routes, such 
as gas-phase con densation or processes on grain surfaces. 



Ceccarelli et al. ( |2002) have also identified carbonate fea- 
tures in the proto-star NGC1 333-IRAS4, and t hey t oo ar- 
gue for a non-aqueous origin. Chiavassa et al. ( 2005 ) have 
identified carbonate features in a range of pre-main se- 
quence environments, from "Class 0" objects to Herbig 
AeBe stars, and identify a need for experimental evidence 
of non-aqueous formation mechanisms. These observations 
demonstrate that carbonates are present in both pre-main 
sequence and highly evolved environments. 

A number of non-aqueous formation routes have been 
proposed to explain the presence of carbonates in circum- 
stellar environments. These include: the processing of CO2- 
rich ices through catalytic react ions on hydrated silicate 



surfaces (Ceccarelli et al. 2002), direct condensation at 



high temperatures from gas rich in CaO and CO2 (Toppani 



2005), and sol id-gas interaction of silicate gr ains and hot, 



gaseo us CO 2 (Kemper et al. [2002 ) . However, Ferrarotti & 
Gail (2005) calculated that direct condensation of carbon- 



ate features in the planetary nebulae (PNe) NGC6302 and 



ate grains would not be viable in such environments due 
to the low temperature (<800K) required for carbonates to 
precipitate from the gas phase. At distances from the cen- 
tral star implied by these temperatures the outflowing gas 
would be significantly diluted, suppressing the formation 
of less stable dust species such as carbonates. In pre-main 
sequence environments carbonates could form through pro- 
cessing on larger planetesimals within the proto-planetary 
disk. Collisions of such bodies would provide a source for the 
carbonate dust observed in these systems. |Rietmeijer et al.| 
( 2008[) support the carb onate formation method suggested 
by i^emper et ai. ( 2002 1 . They have shown experimentally 
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that solid, amorphous CaSiOa, formed though vapour phase 
condensation, as a deep metastable eutectic composition 
which could easily crystal lise as a refractory min eral in as- 
trophysical environments. Rietmeijer et al. ( 2008[ ) proposed 
that these amorphous, Ca-rich silicates could then react 
with gaseous CO2 to produce calcium carbonate (CaCOs) 
via the simple, single-step reaction: 



CaSiOa + CO2 ^ CaCOg + SiOa. 



(1) 



Additionally, Chiavassa et al. ( 2005 ) emphasise the fact 
that both PNe and young stellar objects are known sources 
of X-Ray emission and suggest that carbonates could form 
in both environments through the interaction of X-Rays 
with water-ice on silicate dust grains, increasing the mobil- 
ity of H2O close to that of liquid water. 

1.1.2. Planetary Environments 

Carbonates have lon g been know n to be present within 
chondritic meteorites ( Rubiii|1997 ) and have typically been 
used to infer the presence of liquid water. In the case of me- 
teorites, the presence of carbonate phases is believed to be 
the result of aqueo us alteration of silicate miner als on a 
larger parent body ( Fredriksson fc Kerridge|[l988 ) . A vari- 
ety of carbonates are known to be present on larger bodies 
in the solar system, primarily M ars and Venus ( .Ehlmann] 



distinctive spherical habitat, and although no longer of va- 
teritic composition, they were likely formed inorganically 
by more stable carbonate phases pseudo-morphing an ear- 
lier vaterite assemblage. 

1.2. Vaterite 

The three predominant anhydrous calcium carbonate min- 
eral phases are calcite, aragonite and vaterite. Of these, 
calcite is the most stable, with aragonite and vaterite be- 
ing metastable phases, existing only under a narrow range 
of temperature and pressure conditions. Calcite is therefore 
the most abundant and as such, it is generally calcite that is 
attributed to carbo nate features found within astrophysical 



environments (e.g. Ceccarelli et al. 2002| Chiavassa et al. 



[2005) . However, vaterite has been found in a selection of 



enstatite achondrite meteorites (DuFre sne fc Anders||1962 



Okada et al.||1981| |Vecht fc Ireland||2000 ) 

Naturally forming vaterite is rarely found on Earth, 
predominantly oc curring in association with organic tissue 
( Sutor fc Wooley||1969J and bacteriological biominerahsa- 
tion products (e.g. [Rodriguez-Navarro et al.||2003| Sachez 



IMoral et al.'2003 HFalini et al.|2005|ICacchio et al.|2004[|De 
Muynck ct al. 2008[I Chen et al.|2009(|Zamarreno et al.|2009D 
and freshwater rnollusc shells and pearls (e.g . Wchrmeisterl 

iti et al.|2008]. Naturally 



et al.|2007[pacob et al.|2008,,Soldati et al.|2008) 



et al. 2008 Viviano et al. 2012); and in these cases their occurmg vaterite has been discovered in small amounts 



Kolodny fc Gross|1974[ ), drilhng mud (Friedman et al.|1993 ) 



and in stagnant, natural wat ers, forming as a result of 



spontaneous precipitation (e.g. Rowlands fc Webster||197T 
Grasby|[2003i ) 



presence is also often associated with aqueous alteration within zones of contac t metamorphism ( McConnell| 1960 
during the planets' history. Mars and Venus are both known 
to have abundant CO2 within their atmospheres, albeit at 
vastly different atmospheric pressures, and it has been sug- 
gested that chemical weathering coul d lead to the form ation 
of carbonates through reaction ([T]) ( [Visscher 2009). This 
solid-gas reaction of Ca-rich silicates with CO2 has been 
suggested as a possible method for re gulating the atmo- 
spheric CO2 on Venus (Visscher 2009), with calcium sili- 
cates on its surface reacting with the high surface C O2 sur 
face p ressure of 90 bar 
et al. 11 1968); however 



at a temperature of 750 K (Wood 
this is currently a matter of debate 



Vaterite can be inorganically synthesized in the labo- 
ratory through a number of processes, including precipi- 
tation of a calcium carbonate gel from concent rated cal- 
cium and carbonate solutions ( Andreassen||2005 ) , through 
the decomposition of the hydrated ca,rbonate phase 



ikaite (CaC03.6H20; Tang et al. 2009) or stabilised via 



( [Treiman fc Bullock|2012| and assumes that the abundance 
of Ca-rich silicates exposed on the surface of Venus is suffi- 
ciently high. It has also been argued that this would not be 
a viable buffer of CO2 and there is doubt about the abun- 
dance of woUastonite (CaSi Oa) on the surface of Venu s to 
allow this reaction to occur ( [Treiman fc Bullock||2012[ ). 

Carbonate features have also been reported on the 
surface of Mars using the Compact Reconnaissance 
Imaging Spectrometer for Mars (CRISM) on the Mars 
Reco nnaissance Orbiter ( Ehlmann et al._2008| Viviano et al. 
2012). These carbonates have been found in localised areas 
on Mars, with little evidence for c arbonate minera l s else- 
where on the surface of the planet, van Berk et al.| (|2012[) 
suggests that in Mars' history the partial pressure of CU2 
was much higher than it is today, and that the carbonates 
could have been formed through exposure to CO2 at a par- 
tial pressure of 1 bar over a long period of time (>10^ — 10^ 
years) . 

The Martian meteorite ALII884001 was found to con- 
tain spherical carbonate structures that, due to their well 
defined shape, were initially propo sed as evidence for past 



biomimetic-ba sed precipitation in the pre sence of organic 



additives (e.g. Thompson et al. 2011a|b ), although these 
processes are often involved and, though reproducible, are 
often sensitive to variations in the preparation conditions. 
Once formed, vaterite is metastable and very soluble, more 
so than the calcite and aragonite polymorphs, conve rting 



into calcite in less than 25 hours at room temperature ( Silk 



1970 ) when in contact with water. In the absence of water 
it has been found to transform into calcite at temperatures 



between 730 K and 840 K (Subba Rao 1973 Peric et al. 



The crystal structure of vaterite has been studied in de- 
tail over the last few decades; however the rarity of vaterite 
in nature and its reluctance to form large single crystals 
has lead to a number of uncertainties regarding its crystal- 
lographic structure (see Table [I] for a summary). The main 
report s of t he vate r ite crystal s tr ucture c ome from|Kamhi| 



(1963| and |Meyer| (p60 

gests 



30 th a h exagonal and orthor 



|I969| . |Meyer| ( p60| |1969l ) sug 

d orthorhombic unit cell, wherea! 



life on Mars ( McKay et al.| 1996 ). However, an a l terna- 
tive interpretation proposed by Vecht fc Ireland (2000) 



was that the spherical features were indicative of vaterite 
(metastable CaCOs, see Section 1.2) known to grow in a 



ereas 

Kamhi ( 1963 ) reports a hexagonal unit cell with a promi- 
nent hexagonal pseudo-cell. These cells can be linked, and 
the orthorhombic cell of Meyer can be related to the hexag- 
onal Kahmi cell v ia a transfor mation of lattice parameters 
(see d i scussio n by Tang et aL|[2009) . More recently, ,LeBaiI| 
et al. (2011), in an attempt tont single-crystal data of 
vaterite, propose an Ama 2 space group using a model 
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Table 1. Comparison of proposed crystal structures for 
vaterite from the literature. 



Reference 



Unit Cell Parameters (A) Space Grp. 



Demichelis (2012) 
Mugnaioli (2012) 
Wang (2009) 
Le Bail (2009) 
Kahmi (1963) 
Kamhi (1963) 
Meyer (1960) 



a = 7.12 
a =12.17 
a = 7.29 
a = 8.47 
a = 7.16 
a'= 4.13 
a = 4.13 



b = 7.12 
b = 7.29 
b = 7.16 



7.15 



c = 25.32 P3221 

c = 9.47 C2/c 

c = 25.30 P6522 

c = 4.12 Ama2 

c = 16.98 P&z/mmc 

c' = 8.49 P&z/mmc 

c = 8.48 Pbnm 



based on crystal microtwinning (3 orthorhombic domains 
rotated by 1 20° to produce ps e udohe xagonal symmetry). 



Additionally, Mugnaioli et al. (2012) have used electron 
diffraction techniques, allowing the collection of diffrac- 
tion data from single nanocrystals, to propose a monoclinic 
structure, with a C2/c space group. However, carbonate an- 
ions in vaterite exhibit disorder such that their locations are 



still a matter of debate ( 


Medeiros et al.|2007l|Gebauer et al. 


2009 Jacob et al.||2009 


Wehrmeister et al.||2010D. Density 


Functional Theory (DF 


. ^) calculation by 'Demichelis et al. 



(2012) suggest the previous Pbnm and P6522 structures 
could in fact represent unstable transition states leading to- 
wards a more stable structure with P3221 symmetry. The 
DFT calculations suggest there are at least three distinct 
minima whose energies and activation barriers for intercon- 
version are all within t he accessible range of ther mal energy 
at room temperature. Wang & Becker (2009b), however, 
suggest the P6522 space group for a stable, fully ordered 
vaterite structure, obtained using first principles calcula- 
tions and molecular dynamic simulations. Therefore, the 
uncertainty surrounding the structure of vaterite observed 
in different experiments could be due to differing combina- 
tions of structures involving different symmetries. 

In this paper we report on the formation of vaterite, 
through a solid-gas reaction of amorphous calcium silicate 
powders (CaSiOa) and gaseous CO2 at elevated pressure 
(6-40 bar) and temperature (273 - 1223 K), with the in- 
tention of demonstrating how carbonate phases might form 
in non- aqueous astrophysical environments. We focus on 
carbonation of pure Ca-silicates especially, as Mg-rich car- 
bonates do not fit the observed 92;n m feature, indicative of 
carbonates, in astronomical spectra ( Kemper et al.|[2002 ). 



2. Experimental 

An amorphous calcium-rich silicat e was produced via a so l 
gel method (described in detail by Thompson et al.|2012a ) , 
in which the metal salts CaCl2 and Na2Si03 were combined 
in O.IM solutions to form samples having the stoichiomet- 
ric composition CaSiOa . In order to prevent carbonation of 
the final samples through the interaction with atmospheric 
CO2 , they were stored in gel form under demineralised wa- 
ter until required, at which point the gels were dried in a 
Carbolite HVT vacuum furnace (P <10~''^ mbar), initially 
purged with Nitrogen. The gels were dried at a temper- 
ature of 323 K for approximately 1 hour, producing very 
fine-grained powders (e.g. Figure [T]) with an average grain 
size of <10/im. Upon removal from the furnace samples 
were immediately stored in sealed glass vials under Argon, 
prior to being loaded into the experimental cell. 




TM-1000 0825 



2011/05/13 13:44 



200 urn 



Fig. 1. Scanning electron microscope image of the final 
powder product with composition CaSiOa, vacuum dried 
at 323 K. 



In-situ synchrotron X-ray powder diffracti on (SXPD) 
measureme nts were taken on Beamline 111 (Thompson 
et al.||2009[ ) at the Diamond Light Source, UK using a po- 
sition sensitive detector (PSD) allowing for fast, real-time 



data collection (Thompson et al. 2011a). The 111 beam- 



line receives low divergence, high intensity X-rays from a 
22-pole undulator source, monochromatic X-rays are se- 
lected by a double bounce Si(lll) monochromator, coupled 
with double bounce harmonic rejection mirrors. The PSD 
comprises 18 pixellated Si strip detection modules tiled to 
give a 90° measurement arc. The X-ray wavelength, cali- 
brated against a NIST SRM 640c standard Si powder, was 
0.825590A. 

A high pressure gas cell (Thompson et al. 



2012bl see 
hie 



Figure pi) was used in order to expose the sample to high 
purity (99.98%) CO2 at varying pressures (1-40 bar), and 
heating of the sample was performed using a Cyberstar 
hot air blower (Ramp rate: 10°/min up to 1273 K) located 
beneath the sample capillary. Powder samples were loaded 
into 0.7mm diameter quartz capillaries (maximum pressure 
50 bar) or 0.79mm sapphire tubes (maximum pressure 100 
bar) and held in place with a quartz wool plug. It should 
be noted that, while the sapphire tubes are straight and 
have a constant diameter along their length, the quartz 
capillaries are tapered towards to the (closed) end of the 
capillary; consequently they are usually smaller than their 
nominal diameter. It was therefore found necessary to apply 
greater pressure to the powder to fill the quartz capillaries, 
resulting in tighter packing; as discussed below this has 
implications for the gas pressure required to initiate the 
reaction. The samples were mounted into the high pressure 
cell, aligned on the diffractometer and connected to the gas 
control system ( Parker et al.|2011 ) . Prior to connecting the 
sample to the system the gas supply line was pumped down 
to vacuum to remove air from the system and CO2 was then 
injected to 1 bar. 

Once the sample was mounted the CO2 pressure was 
increased to 6 bar, to ensure the adequate gas penetra- 
tion through the length of the sample to position of the 
X-ray beam. However, due to the slightly smaller diameter 
of the quartz capillary and the tighter packing of the sam- 
ple within the capillary, it was determined that an initial 
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Fig. 2. High pressure gas cell with sapphire tube. X-ray 
beam is at a right angle to the capillary. 



pressure of 20 bar was required in order to obtain adequate 
gas penetration through the sample and observe a reaction 
at the beam position. 

Due to the tighter packing of the sample in the quartz 
capillaries, the gas is prevented from penetrating through 
the complete length of the capillary, instead reacting only 
with the material at the entrance to the capillary, outside 
of the section covered by the beam. In this case, the addi- 
tional pressure required to observe a reaction when using 
the quartz capillary is not believed to be due to any dif- 
ference in the sample but merely represents the additional 
pressure required to all ow the gas to diffuse throu gh the full 



length of the capillary. Thompson et al. (2012a) reported 



the carbonation of the same amorphous samples in a vac 
uum desiccator using ammonium carbonate as the source 
of CO2, supporting the fact that the need for higher pres- 
sure in this study is solely due to the experimental setup, 
and is not necessary for such a reaction to occur. 

The sample was heated steadily from room tempera- 
ture (RT; 298 K) to 1223 K and data were collected at 
5 K intervals. To provide sufficient powder averaging dur- 
ing data collection, the sample capillary was rocked ±15° 
about its length. Every scan with the PSD took 10 seconds. 
For subsequent samples the pressure was increased at am- 
bient temperature, in steps of 10 bar, up to a maximum 
pressure of 40 bar. 



3. Results 

Upon initial injection of CO2, at a pressure of 6 bar, 
into the larger diameter sapphire tube, the formation of a 
crystalline carbonate phase was observed almost instanta- 
neously. Figure [3] shows a comparison of the initial SXPD 
pattern of the CaSiOs sample under vacuum at ambient 
temperature and the first SXPD scan taken as the CO2 
was opened to the sample. Injection of CO2 triggered the 
formation of a crystalline phase. 

Once crystallised the sample exhibits multiple strong, 
well-resolved features, suggesting a well ordered crystalline 
structure. The distinctive groups of features at 10° - 14.5° 
20 , 22° - 24° 29 and 25° - 27° 20 are attributed to vaterite, 
while the remaining features at 12.2°, 15.8°, 19.1° and 20.9° 



29 are due to calcite. Peak search-match software MatchE] 
suggested that the pattern could be fit by a mixture of the 
calcium carbonate phases vaterite and calcite. 

The pressure was held at 6 bar at RT for approximately 

30 minutes. During this time the structure continued to 
evolve, with the peaks steadily growing in intensity; this 
was accompanied by a slow but steady decrease in the CO2 
pressure at the sample. We should note however, that we 
cannot rule out the presence of an amorphous silicate com- 
ponent remaining in the sample; evidence suggests that the 
amorphous silicate component anneals at higher tempera- 
tures, or after prolon ged exposure (^5 hours) at T >1220 K 
(iFabian et al.||2000|). 



3.1. Thermal Evolution 

Heating of the sample was initiated after 30 minutes at 
a pressure of 6 bar. During this time the CO2 pressure 
dropped steadily, an indication of the reaction still occur- 
ring, and was topped up to 6 bar when the pressure dropped 
below 4 bar (after approximately 2.5 hours). The experi- 
ment was then left to run overnight, with scans being taken 
at 5 K intervals (approximately every 3 mins) between RT 
and 1223 K, during which the pressure had dropped and 
settled at 2.6 bar. The results of this are plotted in Figure|4] 
presenting a selection of patterns at roughly 50 K incre- 
ments. These show that, following the initial formation of 
the carbonate phases, they remain relatively stable during 
heating until the temperature exceeds 753 K, at which point 
the carbonate phases start to break down and the sample 
begins to anneal, finally crystallising, at T >1083 K, to 
woUa stonite (CaSiOa). This is discussed in greater detail in 
4.1 



3.2. Pressure Dependence 

Due to the strong features in the diffraction patterns caused 
by the use of a sapphire tube it was decided to perform fur- 
ther experiments using quartz-glass capillaries. As before, 
the sample of amorphous CaSiOa was initially exposed to a 
CO2 pressure of 6 bar; however, scans indicated that no ma- 
jor structural changes occurred upon exposure to the gas. 
The pressure was held at 6 bar for 10 minutes, during this 
time two very sharp, intense peaks appeared at 16.72° and 
33.32° 20 . These correspond to Bragg reflections of calcite, 
however the intensities are much higher than expected, in- 
dicating insufficient powder averaging and therefore likely 
result from the localised carbonation of a small number of 
isolated particles, some of which happen to lie in the Bragg 
condition. As further crystallisation occurs, these particles 
are taken out of the bragg condition and full crystallisation 
occurs at higher pressures. 

The pressure was then increased to 10 bar and held 
for 15 minutes, during which the features at 16.72° and 
33.32° 29 were observed to diminish with time. Before they 
disappeared completely, the pressure was increased to 20 
bar. Initially, there was no observable effect on the struc- 
ture; however, after approximately 15 minutes at this pres- 
sure the steady growth of multiple crystalline features, due 



Crystal Impact http://www.crystalimpact.com/match/ 
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Fig. 3. A comparison of the diffraction patterns for the powdered sample of CaSiOa before (lower curve) and after (upper 
curve) the injection of CO2. The strong feature at 19.97° 29 is due to the sapphire tube. Patterns are offset in the y-axis. 
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Fig. 4. Thermal evolution of CaSiOs exposed to CO2 at 6 bar whilst being heated from room temperature (bottom scan) 
to 1173 K (top scan). Patterns have been offset on the y-axis for clarity and temperature is increasing in steps of 50 K 
from bottom to top. The strong feature at 19.97 29 is due to the sapphire tube. 



to calcite and vaterite phases, was observed. A further in- 
crease in pressure up to 30 bar appeared to accelerate the 
rate of carbonation, increasing the peak intensities, but did 
not produce any further changes to the phase structure. 
Figure [5] shows a selection of SXPD patterns as a function 
of increasing pressure, showing the evolution of the sample. 



4. Analysis 

4.1. Thermal Evolution of CaSi03 



Figure 6(a) shows the result of a Rietveld structure re- 
finement using the TOPAS- Academic software package 



(Coelho 2007) and published lattice parameters for va- 
terite and calcite as starting values. The published lat- 
tice parameters used in the refinement were obtained from 
the Inorganic Crystal Structure Database (ICSD) reference 
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Fig. 5. Evolution of CaSiOs in quartz capillary with increasing pressure of CO2. Patterns have been offset on the y-axis 
for clarity and high backround is due to the quartz capillary. Pressure increases upwards from 6 bar to a maximum of 
30 bar. 



databas^ and are shown in Table [2] The vaterit e phas e 
was fit with the hexagonal pseudo-cell of Kamhi (1963), 
with unit cell parameters of a' = 4.13 A and c' = 8.49 A 
and space group PG^/mmc. The reduced hexagonal cell 
st ructure was ch osen over the larger hexagonal structure 
of Kamhi ( |1963 see Table [l]) as it provided the best fit to 
our data. The inodelled fit is shown in Figure |6(a) along 
with the difference pattern between the overall fit and the 
measured data. The strong peak at 19.97° 20 is due to the 
sapphire capillary and has therefore been excluded from the 
fit. The two phase fit produced agreement factors of R^p 
= 1.78 and Rexp = 1-05 with a goodness of fit (GoF) of 
1.69; the very low R-factors relate to the high background 
present in the data. 

The refined lattice parameters and the calculated weight 
percentage of each phase are shown alongside the published 
values in Table [2j The refined parameters do not differ sig- 
nificantly from the reference values. 

Figure [7] is a plot of weight percentage against tem- 
perature for the vaterite, calcite and woUastonite phases, 
with values taken from a batch refinement of selected data 
sets at approximately 20 K temperature intervals. This pro- 
vides a more detailed view of what is happening, with re- 
gards to composition of the sample at higher temperatures 
(>800 K). From this we can clearly see that, at lower tem- 
peratures, the vaterite and calcite phases remain fairly sta- 
ble, with the relative weight percentage of vaterite increas- 
ing by just 3% over a period of approximately 5 hours. 
However, once the temperature reaches 723 K the percent- 
age of vaterite begins to fall, being replaced by an increas- 
ing amount of calcite. The weight percentage of vaterite 
continues to fall gradually until calcite supersedes it at a 
temperature of 770 K. At a temperature of 840 K, woUas- 



Table 2. Published and refined lattice parameters for va- 
terite and calcite phases. Published values obtained from 
the ICSD database. Error values stated are for the final fig- 
ures. Values given are from samples at room temperature 
and initial formation pressure; 6 bar for sapphire, 20 bar 
for quartz. 



Vaterite - 


Space Group: PGa/mmc 






Capillary 




a' (A) 


c' (A) 


Wt % 




Published 


4.13 


8.49 




Sapphire 


Refined 


4.1226 (1) 


8.4653 (3) 


65.5(4) 


Quartz 


Refined 


4.1227 (4) 


8.462 (1) 


50.2(4) 


Calcite - 


Space Group: 


R-3c 






Capillary 




a (A) 


c(A) 


Wt % 




Published 


5.05 


17.32 




Sapphire 


Refined 


4.9869 (1) 


17.0492 (4) 


34.4(4) 


Quartz 


Refined 


4.9844 (4) 


17.058 (1) 


49.7(4) 



^ F IZ Karlsruhe jhttp: / /cds. dl.ac.uk/cds/datasets/crys/icsd/ 
|llicsd.html| 



tonite begins to crystallise and can be included in the fit, 
this phase then quickly dominates the sample as it contin- 
ues to crystallise at higher temperatures. Figure [4] indicates 
that at a temperature of 973 K, the carbonate peaks have 
almost disappeared, appearing as weak, broad features, in- 
dicating that the carbonate phases are not stable at such 
temperatures and have broken down, this is also evidenced 
in Figure [7] as the amount of vaterite and calcite within 
the sample at this stage are less than 10 %. This leads to 
the sample fully annealing at a temperature of 1103 K to 
fully crystalline CaSiOa (woUastonite). Alternatively, plot- 
ting lattice parameter a against temperature for the va- 
terite phase indicates a slight thermal expansion between 
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Fig. 7. Weight percentage against temperature for calcite 
and vaterite phases of initial CaSiOa sample in a sapphire 
capillary, at a pressure of 6 bar. Where error bars are not 
visible they are smaller than the plotted data points. 

RT and 943 K, this provides a thermal expansion coefficient 
of-SxlO^^K-i (see Figure[8|. 

4.2. CaSiOs in Quartz Capillary 

The crystalline phase formed at a pressure of 20 bar on the 
second run was identified as also being a combination of cal- 
cite and vaterite. Refinement of the data using initial crys- 
tallographic parameters obtained from the ICSD data base, 
confirmed this. A typical fit is shown in Figure 6(b) with 
the initial and refined parameters and calculated weight 
percentages listed in Table [2] 




4.05 



300 400 500 600 700 800 900 1000 
Temperature (K) 



Fig. 8. The effect of temperature on lattice parameter a 
for the vaterite phase. Sample in a sapphire capillary at a 
pressure of 6 bar. 

The refined values are in good agreement with published 
values and the refined values for the sapphire capillary sam- 
ple. Figure [9] shows the weight percentages of the two phases 
plotted against pressure and it can be concluded that once 
the two phases have been formed, they remain relatively 
stable. An increase in pressure does not appear to have a 
significant effect on the composition of the two phases. 

5. Discussion 

5.1. Formation of Vaterite 

As discussed in the introduction vaterite is, with respect to 
aragonite and calcite, the least stable of the anhydrous cal- 
cium carbonate phases. Its formation in the present study 
is therefore interesting, either as a possible component of 
dust grains, or as a possible metastable precursor phase 
for the calcite grain component curr ently attributed to the 
92/im feature (Kemper et al. 2002). We have shown that 



100 
90 
80 
70 
60 
£ 50 
40 
30 
20 
10 



"T 1 1 1 1 1 T" 



- Calcite 

- Vaterite 



_l I I I I I I L_ 



20 



25 30 
Pressure (Bar) 
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Fig. 9. Weight percentage against pressure for calcite and 
vaterite phases at room temperature. Initial sample of 
CaSiOa in quartz glass capillary. Where error bars are not 
evident, they are smaller than the plotted data points. 



the calcium carbonate phase vaterite can be formed through 
the solid-gas interaction of amorphous calcium-rich silicates 
and gaseous CO2 at an elevated pressure of 6 bar. While the 
kinetics of the transformation in to calcium carbonate are 
dependent on the initial packing density of the powdered 
CaSiOa sample, as discussed in Section 2, the selection of 
vaterite over the more thermodynamically stable calcite is 
likely to be due to material dependent factors. In the fol- 
lowing paragraphs we draw on the results of experiments 
regarding vaterite formation in other fields to show that 
the preferential stabilisation of vaterite in the present case 
likely arises from the presence of Si02 anions. 

Due to its industrial and commercial potential, there is 
a large body of literature that shows vaterite can be sta- 
bilized (often temporarily) in the presence of various dif- 
ferent ions or organic and inorganic additives under widely 
differing experimental conditions, most of which are con- 
cerned with in vitro formation (e.g. Wray &: Daniels||1957 



Nassrallah-Aboukais et aT] 1998[ |Nothig-Laslo fc Brecvic 



1999; Han et al. 2006': 'Nebel fc Epple ||2008[ |Wang et al." 
,2009; Fernandez-Daz 2010; Gomez-Morales et al.||2010| ) , or 
as a post-formation coating (Vogel et al. 2009p . The ori- 
gin of these stabilisation effects likely lies in the fact that 
vaterite, aragonite and calcite all have very simi lar values 
for their free energy (0.5-3 kJ/mol differences; Plummer 
fc Busenbergl|T982| |Wolf et al.|[l996l [2000; Baitalow et al 



1998). Relatively small changes in either the surface area 



i.e. small particle size), or in the concentration of impu- 
rities, can change their stabilisation properties such that 
the ord er of their therm odynamic stabilities becomes in- 
verted (Navrotsky 2003 |2004) . Similar thermodynamic 
crossovers are observed for several metal oxide systems, in- 
clu ding the Ti02 po lymorphs (rutile, brookite, a nd anatase; 
iRanade et al.|[2002 ) and the iron oxide phases (Navrotsky 
|et al. 2008p and can mean that the least stable phases, 
in the present case vaterite, are often preferentially formed. 
However, many experimenters (see papers cited above) have 
reported that, although outcomes are highly reproducible, 
the actual outcome is also highly depende nt on changes in 



conditions in a non-systematic way (e.g. |Pai et aL 2009 



Ren et al. 2009), suggesting that kinetic factors may also 
be important. It is also relevant to note Ostwald's Law of 
Phases, which states that, due to the lower energy barrier 
of more disordered states, the pathways to final crystalli- 
sation of the most stable phase will pass through all less 
stable phases in order of increasing stability. 
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For precipitating carbonates a hydrated amorphous cal- 
cium carbonate phase (ACC) acts as a precursor, which 
can then crystallise to a more stable form. Carbonate for- 
mation via this route is believed to follow: ACC — » va- 
terite — >■ aragonite — »■ ca lcite (e.g. Rodriguez-Bianco et al. 

al. 2012 and references therein). In this 
Telieved that due to the local coordina- 
of Ca in the nanoframework structure 



2011 Bots et 



instance, it is 
tion environment 
of ACC being similar to t hose in the crystalline CaCOa 



phases (Becker et al. 2003), the dehydration and conden- 



sation 01 the ACC structure results in the formation of va- 



terite without major structural rearrangements ( Gebauer 
erar][20T0l [Goodwin et~ar][20T0l [Rodriguez-Bianco et a 



20111. In the case of carbonate formation by biomineral- 
isation, th i s ofte n occurs in the presence of silicates and 
Gal et al. (20101 have suggested that a mechanism of ge- 
ometric frustration may operate, whereby the presence of 
the tetrahedral silicate ion in the flat carbonate lattice hin- 
ders the organi sation of the carbonate into crystalline poly- 
morphs; while Lakshminarayanan & Valiyaveettil (20031 
previously observed a signiflcant expansion of the calcite 
lattice with added silicate anion concentration, attribut- 
ing this to the incorporation of monomeric, oligomeric, and 
polymeric forms of silicate anion species into the calcite lat- 
tice. Such incorporations would increase the calcite lattice 
energy, making it less stable and therefore less favourable. 

Calcium carbonate formation in the presence of var- 
ious s ilicate s and oxides was also investigated by |Lin[ 
( 2009 1 , who identified a correlation between the sur- 



et al. 



face charge provided by the silicate/oxide and the choice 
of phase stabilised. Negative surface charge favoured the 
metastable vaterite and aragonite phases, while positive 
charges produced calcite. In crystalline silicates, the pos- 
sible charge sites are limited by the crystal structure: 
planar regions which have a negative structural charge 
and edge regions whose charge arises from the protona- 
tion/deprotonation of surface hydro xyls at the side of stack- 



ing structures (sheets, chains etc.). Lin et al.[ ([2009[) found 
amorphous Si02 did not follow the surface charge rule and 
always stabilised calcite. However they attributed this to 
the gel phase that forms when Si02 is exposed to water, 
suggesting that the hydrated porous nature of the gel is 
not compatible with epitaxial growth and therefore no tcm- 
plating effect is exerted on carbonate nucleation, allowing 
calcite to form rather than vaterite or aragonite. 

Assuming all the Ca atoms eventually participate in 
CaCOa, for carbonate formation to proceed via the gas- 
solid reaction in CaSiOa necessitates the formation of Si02, 
providing a free ionic species, while the absence of liquid 
water will negate the formation of a gel phase, allowing 
Si02 tetrahedra to either be partially inc orporated into the 
the carbonate structure, as suggested b y Gal et al. (2010) 
and Lakshminarayanan & Valiyaveettil ( 2003p ; or, to exert 
a templating su rface charge effec t on the forming carbonate 
as suggested by Lin et al. ( 2009 1 . Since carbonation will oc- 
cur at the Ca sites within the silicate structure, the initial 
small size of the carbonate could also contribute to the car- 
bonate free energy to produce a Navrotsky stability inver- 
sion and similarly favour the stabilisation of vaterite. The 
fact that some calcite is also initially formed which, in the 
second experiment, was observed to be unstable possibly 
points to contributions such as particle size and remaining 
silicate "impurities" also influencing phase stabilisation pri- 
orities. We also acknowledge that kinetic factors may play 



a role in determining vaterite stabilisation and subsequent 
growth, particularly as our initial gas injections were per- 
formed at room temperature and relatively high gas pres- 
sures to overcome sample packing density. However these 
effects are also likely to favour the formation of disordered 
structures, i.e. vaterite rather than aragonite or calcite. 

The vaterite formed in both cases appears to be of 
the same structure, both refine w ell using the hexagonal 
pseudo-cell values of Kamhi ( 1963 ) . Table l2]lists the refined 
lattice parameters for vaterite at room temperature from 
both sources. The values are in good agreement with each 
other, suggesting that, although formed independently, the 
structure of the vaterite is reproducible. 

5.2. Relevance for Cosmic Dust 



A ban d near 90/zm was first discovered by [Kemper et al. 
( [2002[ ) in the spectrum of the PNe NGC 6302, followed 
by detections in the class proto-star IRAS 16293 and 
two HII regions (Ceccarelli et al.[[2002| [Takashi & Yoko 



[2003); a subsequent survey of low mass pr oto-stars showed 



2005). The circumst e llar fe ature at 
By" Ceccarelli et al. 



the fe ature to be common in such regions ( Chiavassa et al 

92/im was attributed 
calcite, who ruled out the 



^ ([20021 to 

metastable aragonite phase on the basis that its 92//m fea- 
ture is too weak and that its formation nominally requires 
high pressures and temperatures. Whether vaterite exhibits 
an observable Far Infrared (FIR) feature, or indeed could 
be a possible candidate for the carrier of the 92/im fea- 
ture, is difficult to assess. There are, to our knowledge, no 
published FIR spect ra for vaterite (e.g revie w of carbon- 



ate spectroscopy by Brusentsova et al. 2010), likely made 



difficult due to a combination of its low thermal stability 
in respect of producing polyethelene pellets and the dif- 
ficulty of producing or obtaining pure specimen. However 
the lattice mode vibrational region is accessible by Raman 
spectroscopy. The Raman spectroscopic signatures of car- 
bonates are found in three wavenu mber regions: 1500-1000 , 
1000-500, and 500- 100 cm"! (e.g. |Scheetz fc White[[T977 



Ed wards et al. 2005 ). In general, bands at frequencies above 
500 cm~^ are due to the internal motions of the molecu- 
lar carbonate ion (internal modes), and those below 500 
cm^^ are due to motions involving the entire unit cell 
(lattice modes). In this region the carbonate polymorphs 
exhibit distinct and characteristic spectra l features (e.g. 
Edwards et al. 2005 Carteret et al.[ 2009 1 . In particular 
vaterite exhibits a strong broad feature at ^ 100 cm^^ 
and is illustrated in Figure 10 which shows Raman spec- 
tra for a number of different vaterite samples produced 
via a bi omimetic stabilisation process (Thompson et al. 
2011a|b I that involved the hydrolysis of urea in the presence 



of amino acid. In these preparations varying quantities of 
vaterite and aragonite were precipitated along with a minor 
quantity of calcite. Each polytype is normally distinguish- 
able by its crystal morphology: needles for aragonite, cubes 
for calcite and spherules for vaterite. For the amino acid, 
leucine, aragonite exhibited as splintered or branched nee- 
dles, calcite as stepped/intersecting cuboids and vaterite as 
fiower-like platelets. By focussing the laser spot on the dif- 
ferent morphological structures, pure-phase Raman spectra 
were obtained (Horiba LamRam 800, 532 nm NdYAG laser, 
600 line grating, 10 x objective). In the figure, spectra are 
shown for five samples precipitated under differing condi- 
tions. As can be seen in Figure [Tol the Raman spectra for 
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Table 3. Factor group analysis a nd selection rules for cal- 
cite (after Decius &: Hexter|[l977 ). 



Dm 




A2g 


Eg 




A2u 


Eu 


Ca 2b 











1 


1 


2 


C 2a 





1 


1 





1 


1 


O 6e 


1 


2 


3 


1 


2 


3 


Nt 


1 


3 


4 


2 


4 


6 


Ta 











1 


1 


T 





1 


1 


1 


1 


2 


R 





1 


1 





1 


1 




1 


1 


2 


1 


1 


2 


Activity 


Raman 




Raman 




IR 


IR 



Notes. Nt = Total number of unit cell modes; Ta = acous- 
tic branch modes; T — optic branch translatory modes; R = 
librational modes; F"*'' = internal modes. 

Table 4. Factor group a nalysis and selection rules fo r va- 
terite lattice modes (after [Anderson &: Brecevi3||1991 1. 





A2u 


Eiu 


Alg 


Elg 


E2g 


-pvib 


1 


4 


2 


2 


4 


R 


1 


2 


1 


4 


2 


Activity 


IR 


IR 


Raman 


Raman 


Raman 



Notes. R — librational modes; F 
modes. 



vib / T^vib 

*--^Ca^ 



CO3 



— internal 




both IR and Raman active (e.g. White 2005 Hernandez 
eraL][T997l [Duley et al.||2005[ [Hopkinson et al.||2008[ |Abe 



fc Shigenari||2011 ). In the case of vaterite, the uncertainty 



regarding its unit cell structure and the level of disor- 
der ma ke factor grou p analysis less certain. However, us- 



ing the Meyer] ( 1969 ) PQ^/mmc space group, Anderson & 
'Brecevic' f 1991') identified the selection rules listed in Table 
[4] which suggests there are eight IR active and 15 Raman 
active vaterite lattice modes. Based on assignments to mea- 
sured bands observed at higher frequencies Anderson and 
Brecevic predict three IR active bands (corresponding to 
librations + 2£'i„) should lie below 200 cm~^. 

Although none of the forgoing is a "proof-positive" argu- 
ment that vaterite possesses a FIR band near 92/^m, it does 
suggest that measurement of vaterite in the FIR should be 
addressed as a matter of priority, particularly as the disor- 
der inherent in the vaterite structure might also provide a 
natural mechanism to account for the observed variations 
in the carbonate feature width and position. 

Although the operational pressures required for the 
present work arise from the need to provide adequate dif- 
fusion of CO2 along the length of the sample capillary, we 
believe similar interactions involving silicates in circumstel- 
lar outfiows could occur at much lower pressures. However, 
even if this is not the case the formation of carbonate via 
this route may still be directly relevant to proto-stellar and 
proto-planetary envionments, where carbonates formed at 
higher pressure on, or in, planetesimal objects could be re- 
leased as dust grains via disruptive coUisional events. 




300 

Raman shift (cm"') 



Fig. 10. Raman spectra for samples of CaCOs vaterite pro- 
duced by biomimetic synthesis, showing a strong lattice 
mode feature at ~ 100 cm~^. Spectra offset in y-axis di- 
rection for clarity and y-axis scale is arbitrary. 



these vaterites show a strong feature in the region of 110 
cm~^('^91/im). 

The lattice-mode assignmen t for calcite is well estab- 
lished ( Decius fc Hexter 1977) with the selection rules 
shown in T'able[Sj T'he coupling pvhich leads to factor-group 
selection rules is relatively weak such that lattice modes 
of a given basic origin (e.g. libratory) will be f ound in 
similar positions in both IR and Raman spectra (Adams 
& Williams 1980). In addition, structural disorder, which 
shows up as broadening of Raman feature widths in the 
lattice mode region, is a symmetry reducing effect which 
can in principle lead to a relaxation of the Raman and IR 
activity selection rules such that certain bands can become 



6. Conclusions 

We have observed the formation of the rare calcium carbon- 
ate polymorph vaterite through the solid-gas interaction of 
amorphous silicates, produced as cosmic dust analogues, 
and gaseous CO2 at pressures >6 bar. In-situ synchrotron 
X-ray powder diffraction was used to observe the structural 
evolution of the material during its formation and subse- 
quent heating to a temperature of 1173 K. We determined 
a range of stability between 6-30 bar and 298 K - 923 K, 
a much wider range than previously seen for vaterite. This 
has relevance to the surface cornpositions of Venus and an- 

2]), where the CO2 pressure 



cient Mars (van Berk et al. 2012^ 



at the surface is sufficiently high for this kind of reaction 
to occur. Vaterite has been regarded as important due to 
the fact that it is rarely formed naturally and therefore 
its presence is usually taken as being indicative of biogenic 
processes. 

Our alternative, inorganic method provides evidence 
that vaterite found in Martian meteorites and on the sur- 
face of Venus and Mars are not necessarily proof of bioac- 
tivity but could have formed naturally through a processes 
similar to those that we describe here. 
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Fig. 6. Sample rietveld refinements of powder patterns exhibiting vaterite and calcite peaks. The red curves show the 
fit, the grey curves show the difference curve and the markers indicate the positions of the individual vaterite (blue) and 
calcite (green) peaks, a) Initial sample exhibiting vaterite features in sapphire capillary. The sapphire capillary exhibits 
a strong peak at 19.97° 26 which has been removed from the plot, b) Powder pattern exhibiting vaterite peaks at a CO2 
pressure of 30 bar in quartz capillary. 
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